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ABSTRACT 

The Sunyaev-Zel’dovich effect (SZE) observable-mass ( Y-M ) scaling relation is a promising technique for 
obtaining mass estimates for large samples of galaxy clusters and holds a key to studying the nature of dark 
matter and dark energy. However, cosmological inference based on SZE cluster surveys is limited by our 
incomplete knowledge of scatter, and evolution in the Y-M relation. In this work, we investigate the effects of 
galaxy cluster mergers on the scaling relation using the Omega500 high-resolution cosmological hydrodynamic 
simulation. We show that the non-thermal pressure associated with merger-induced gas motions contributes 
significantly to the scatter, and evolution of the scaling relation. After the merger, the kinetic energy of merging 
systems is slowly converted into thermal energy through dissipation of turbulent gas motions, which causes the 
thermal SZE signal to increase steadily with time. This post-merger evolution is one of the primary source 
of scatter in the Y-M relation. However, we show that when the missing non-thermal energy is accounted 
for, the resulting relation exhibits little to no redshift evolution and the scatter around the scaling relation is 
~ 20-30% smaller than that of the thermal SZE signal alone. Our work opens up a possibility to further 
improve the current robust mass proxy, Y, by accounting for the missing non-thermal pressure component. We 
discuss future prospect of measuring internal gas motions in galaxy clusters and its implication for cluster-based 
cosmological tests. 

Subject headings: cosmology: theory - galaxies: clusters: general - galaxies: clusters: intracluster medium 


1. INTRODUCTION 

Galaxy clusters are the largest gravitationally collapsed ob¬ 
jects in the universe. Their formation and evolution are crit¬ 
ically sensitive to details of the expansion history of the uni¬ 
verse. Accurate measurements of their abundance as a func¬ 
tion of cluster mass and redshift can therefore provide impor¬ 
tant constraints on cosmological parameters such as Qm, Ha, 
erg and vv’o (see, e.g. Vikhlinin et al. j2009| Mantz et al. |2010Jl. 
However, the use of galaxy clusters as cosmological probes is 
currently limited by our ability to relate observable properties 
and cluster mass through so-called observable-mass relations. 

In recent years, the Sunyaev-Zel’dovich effect (SZE) has 
emerged as a promising avenue for both detecting clusters and 
obtaining robust mass estimates for hundreds of clusters in 
the high-redshift universe (Carlstrom et al. [2002] for review). 
The SZE is a small distortion in the cosmic microwave back¬ 
ground (CMB) spectrum caused by the scattering of CMB 
photons off a distribution of high-energy electrons in dense 
stmetures such as clusters of galaxies (Sunyaev & Zeldovich 
|1970| |1972[ Birkinshaw |1999[ for review). The SZE has the 
unique property that its signal is independent of redshift and 
does not suffer from cosmological surface brightness dim¬ 
ming, in contrast to X-ray observations which have been the 
primary observational tools for studying the hot X-ray emit¬ 
ting intracluster medium (ICM) in galaxy clusters. Moreover, 
the integrated SZE signal, Y, is only weakly dependent on 
nonlinear astrophysical processes (e.g. cooling flows and en¬ 
ergy feedback from supernova and supermassive black holes) 
which play an important role in determining the properties of 
cluster cores and significantly affect their X-ray luminosity, 
making the SZE an excellent probe of cluster mass. 


The current generation of SZE cluster surveys, conducted 
with the Atacama Cosmology Telescope (ACT), the South 
Pole Telescope (SPT), and the Planck space mission, have ob¬ 
tained large (~ 10 4 ) and statistically representative samples of 
massive clusters_across the full sky (Hasselfield et al. 2013] 
Reichardt et al. |2013[ Planck Collaboration et al. |2014b[ ) out 
to a median redshift significantly higher than that of X-ray se¬ 
lected cluster catalogs. However, to realize the full statistical 
power of SZE surveys, systematic uncertainties in the Y — M 
relation must be controlled at a level comparable to statistical 
uncertainties. For example, Planck Collaboration has recently 
reported a tension between the cosmological constraints from 
CMB and galaxy clusters (Planck Collaboration et al. [2014a]), 
possibly indicating an unexpectedly large bias in the calibra¬ 
tion of the Y-M relation (but see also Rozo et al. |2014) >. A 
detailed study of the Y-M relation is, therefore, particularly 
timely and important. 

Theoretically, modern hydrodynamical simulations suggest 
that a SZE signal integrated to a sufficiently large fraction 
of cluster volume is a robust mass proxy with a scatter of 
|2005| Nagai [2006] Stanek et al. |20l0 , 
2). Specifically, while the normalization 
is sensitive to cluster astrophysics such as gas cooling, star 
formation, and energy feedback from stars and active galactic 
nuclei (AGN), the scatter and evolution of the Y-M relation is 
fairly insensitive to the poorly understood cluster astrophysics 
(Nagai 2006[ Battaglia et al. |2012| Fabjan et al. 20111. Cluster 
merger dynamics, on the other hand, has been shown to play 
an important role in determining the scatter and evolution of 
the Y-M relation (Wik et al. 120081 Yang et al.[2010; Krause 
et al. |20l2| ). 

This picture appears to be supported by recent observations 
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which showed that dynamically disturbed clusters have in¬ 
creased scatter around the relation compared to their more 
relaxed counterparts (Sifon et al. |2013j l. To harness the po¬ 
tential of the SZE surveys, we must understand the origin of 
scatter, and evolution in the Y —M relation and devise a tech¬ 
nique to control their systematic uncertainties to better than a 
few percent for a wide range of masses and redshifts. 

In this work, we use the Omega500 high-resolution cosmo¬ 
logical simulation to demonstrate that the scatter and evolu¬ 
tion in the Y—M relation are governed primarily by the evolu¬ 
tion of the non-thermal pressure provided by internal gas mo¬ 
tions generated by major mergers. To characterize the effects 
of mergers and dynamical states of clusters, we ( 1 ) analyze 
non-radiative simulations in which the ICM is heated by dom¬ 
inant gravitational heating processes (such as shock and com- 
pressional heating) alone and ( 2 ) use detailed merger trees to 
track each cluster’s formation history and directly examine 
the effects of major mergers on the scaling relation. We first 
show that large deviations from the scaling relation in unre¬ 
laxed systems are due to a large influx of non-thermal pressure 
from random gas motions generated by major mergers. We 
also find that the scaling relation departs from self-similarity 
at high redshift where mergers are more ubiquitous. By ac¬ 
counting for the missing non-thermal energy contribution, we 
demonstrate explicitly that the thermal plus non-thermal pres¬ 
sure exhibits remarkable regularity and reduced scatter at all 
redshifts, indicating that major mergers are a dominant source 
of uncertainty in this relation. 

This paper is organized as follows: in Section[2]we describe 
the thermal SZE and present relevant scaling relations pre¬ 
dicted by the self-similar model, as well as previous results on 
the evolution of thermal pressure support after major mergers; 
in Section [3] we describe the Omega500 cosmological simu¬ 
lation; in Section[4]we present our analysis and interpretation 
of the evolution of SZE-mass scaling relations; finally, in Sec¬ 
tion [5] we summarize our findings and discuss their implica- 
tionsfor current and future SZE cluster surveys. 

2. THEORETICAL FRAMEWORK 
2.1. Thermal SZE 

The thermal SZE is a distortion in the CMB spectrum 
produced by the inverse Compton scattering of CMB pho¬ 
tons off free electrons in dense structures such as clusters 
of galaxies (Sunyaev & Zeldovich |1970[ |1972[ ). This scat¬ 
tering causes a small change in the mean photon energy as 
Av/v ~ k B T /( m e c 2 ) ~ 10~ 2 , where v is the frequency of the 
CMB photon and m e is the electron rest mass. As a result, 
the intensity of the CMB spectrum increases in the high- 
frequency (Wien) end and decreases in the low-frequency 
(Rayleigh-Jeans) tail. The corresponding brightness fluctu¬ 
ation in the CMB is of order ~ 10 -4 , about an order of mag¬ 
nitude larger than the fluctuations from primary anisotropies. 

The change in the CMB specific intensity at a frequency 
v caused by the thermal SZE is proportional to the line-of- 
sight integral of the number density (n e ) and temperature (T e ) 
of electrons, and is given by A/„/7 cmb = fu(x)g,Ax)y, where 
x = hv/k B T C uR, fJx) = x(e x + l)/(e* - 1 ) -4, and g v (x) = 
x A e x j(e x — l) 2 . The dimensionless Comptonization parameter 
y is defined as 

y = k -^L [n e {l)T e {l)dl (1) 

m e c z J 

where c is the speed of light, and ctj is the Thomson cross¬ 


section. This equation corresponds to the integrated line- 
of-sight thermal pressure of the intracluster gas, which is 
approximated as P = n e k B T in the ideal gas limit. The 
corresponding change in the CMB temperature is given by 
A7j,/7cMB = fv(x)y. In the Rayleigh-Jeans limit (u <C 
200GHz), ATjT CM B=-2y and A I v = Qk B v 2 /c 2 )AT„. 

Let us now consider the SZE signal arising from a clus¬ 
ter located at redshift z. The SZE flux integrated across the 
surface of a cluster is defined as the integrated Compton-y pa¬ 
rameter Ysz- 

Y= f ydSl=-±-( [ n e {l)T e (l)dV , (2) 

Jn dy(z.) \ m e c z J J v 

and dV. = dA/d\{z ) is the solid angle of the cluster subtended 
on the sky (the integral is taken over the volume of the clus¬ 
ter), aU(z) is the angular diameter distance to the cluster, dA is 
the area of the cluster on the sky, and dV is the cluster volume. 
Y measures the total thermal energy of a cluster. 

The thermal SZE signal is linearly sensitive to gas mass 
M gas = f gas M and mass-weighted temperature T m : 

Y oc fgpsMT m . (3) 

where f gas is the gas mass fraction and M is the total cluster 
mass. In this work, we use a spherically integrated Y within 
a radius of interest around the cluster center, rather than the 
projected Y derived from observations. This enables our work 
to focus on the direct effects of mergers and dynamical state 
on Y. 

2 .2. Self-similar Scaling Relations 

In the absence of cooling and heating processes, clusters 
are expected to scale self-similarly (Kaiser ] 198~6| . The self¬ 
similar model predicts that the temperature of the gas scales 
with the cluster mass as 

MtxT% 2 E-\z) (4) 

where M = 4Trr 3 A A c p crit /3 is the halo mass enclosed within 
rA, defined as a radius of spherical volume within which the 
mean density is A f times the critical density, p CI j t , at that red¬ 
shift (Bryan & Norman [l998| . E(z) is the redshift-dependent 
Hubble parameter, defined as H(z ) = 100/z£’(z) km s _1 Mpc , 
and it is given by E 2 (z) = H m (1+z ) 3 + Ha for a flat cosmology. 

Inserting Eq. [4] into Eq. [3] we obtain the SZE scaling rela¬ 
tion predicted by the self-similar model, 

Y cx f glis M 5,3 E 2 ^(z), (5) 

where f gas is taken to be constant in the self-similar model, 
such that Y cx M 5 ^ 3 E 2 ^(z). In practice, we fit a straight line 
to the InT-lnM relation by minimizing \ 2 , where the best-fit 
relation is given by 

lnY fit = lnA u +a\n ^ 1 Q i 4 ^i M ) + /? ln E(z), ( 6 ) 

where In A 14 is the normalization constant of the best-fit rela¬ 

tion at M = lO 14 /r'M 0 , and a is the best-fit slope of the Y-M 
relation. Note that the self-similar model predicts a = 5/3. /? 
is the redshift dependence of the relation and is fixed to the 
self-similar value of 2/3. 

Throughout this paper, we consider Y, mass and other clus¬ 
ter properties within two commonly used radii, defined by the 
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total matter overdensity they enclose and chosen for their rel¬ 
evance to X-ray and SZ observations. We use radii rsoo (: and 
f 200 c, enclosing overdensities of A e = 500 and 200 times the 
critical density of the universe, p cm and radius r 2 oom enclos¬ 
ing overdensities of A m = 200 times the mean density of the 
universe. 


2.3. Evolution of Thermal and Non-thermal Pressure Content 

The evolution of thermal and non-thermal pressure support 
after mergers has been studied extensively by Nelson et al. 
( |2012| ). Nelson et al. ( |2012| hereafter N12) have shown that 
during a major merger, non-thermal pressure from random gas 
motions provide about a third of the total gas pressure in the 
ICM. As clusters relax, the non-thermal pressure support falls 
off rapidly, ultimately comprising only a residual 10% con¬ 
tribution. This decay of non-thermal pressure support is due 
to the conversion of random velocities to thermal energy via 
shocks and dissipation. As a result, the total energy content 
(thermal and non-thermal) of the ICM remains nearly constant 
after a merger - a decrease in random kinetic energy corre¬ 


sponds to an increase in thermal energy (Vazza et al. 2011 
Nelson et al. |2012j ). In this paper, we consider properties of 
the ICM measured under the assumption of spheric al sym - 
metry following the methods described in Lau et al. (|2009|>. 
The total energy content of cluster gas has two components: 
the thermal energy from the hot ICM and non-thermal energy 
from gas motions, herein defined respectively as 


eth = J Pth dV (7) 

Em = j ^Pgas(c , 2 + of) dV (8) 

where P,i, is the thermal pressure, p gas is the gas density, and 
a ]: and of are radial and tangential dispersions measured in 
spherical shells, using the peculiar velocity of the cluster dark 
matter within r S oo r to define a rest frame. However, we have 
verified that our results are insensitive to the choice of rest 
frame. Note that a complete description of non-thermal en¬ 
ergy includes rotation, streaming and cross terms. However, 
these contributions are small compared to that of random mo¬ 
tions (Lau et al. |2013| > and neglected in this work. 

As stated above, the Y—M scaling relation is dependent on 
a tight correlation between the thermal properties of the ICM 
and the mass of the cluster. However, as shown in N12, a 
large amount of non-thermal pressure is injected into a clus¬ 
ter during a major merger which only thermalizes over a few 
Gyr timescale, resulting in a prolonged deviation from this 
relation. To characterize the contribution from the missing 
non-thermal energy contribution, we calculate a combined Y, 

Y th +Y nt = S -±^Y th , (9) 

Sth 

where K,/, and Y nt are the thermal and non-thermal components 
of Y . Y,i, + Y„ t is proportional to the total energy content of the 
ICM and therefore more closely correlates with the total mass 
of the system. Note that spherical symmetry for the gravita¬ 
tional potential and steady state are assumed in deriving the 
expression for non-thermal energy, and all the physical quan¬ 
tities at a given radius are averages over a radial shell (Lau 
et al. 12013) . 


3. SIMULATIONS 

3.1. Hydrodynamic Simulations of Galaxy Clusters 

We analyze the Omega500 simulation of galaxy clusters 
presented previously in Nelson et al. ( 2014a[ >. We refer the 
reader to that paper for more details. We briefly summarize 
the relevant parameters below. 

In this work, we analyze the Omega500 high-resolution 
cosmological simulation of massive galaxy clusters in a flat 
ACDM model with WMAP five-year ( WMAP5 ) cosmologi¬ 
cal parameters: Hm = 1 - Cl a = 0.27, f lb = 0.0469, h = 0.7 
and < 7 g = 0.82, where the Hubble constant is defined as 
100/z km s -1 Mpc 1 and erg is the mass variance within spheres 
of radius 8 /z _1 Mpc. The simulation is performed using 
the Adaptive Refinement Tree (ART) ;V-bodv+gas-dynamics 
code (Kravtsov |1999[ Kravtsov et al. |2002[ Rudd et al. |2008| ), 
which is an Eulerian code that uses adaptive refinement in 
space and time, and non-adaptive refinement in mass (Klypin 
et al. | 2001 j > to achieve the dynamic ranges to resolve the 
cores of halos formed in self-consistent cosmological simu¬ 
lations. The simulation volume has a comoving box length of 
500/z -1 Mpc, resolved using a uniform 512 3 grid and 8 levels 
of mesh refinement, implying a maximum comoving spatial 
resolution of3.8/z _l kpc. 

We selected clusters with M 5 oo c > 2.9 x 1O 14 /z-'M 0 and 
performed a simulation where only the regions surrounding 
the selected clusters are resolved. The resulting simulation 
has effective mass resolution of 2048 3 surrounding the se¬ 
lected clusters, allowing a corresponding mass resolution of 
1.09 x 10 9 /z _1 Mg. The current simulation only models grav¬ 
itational physics and non-radiative hydrodynamics. Nagai 
(2006) has studied the effect of cooling and star formation on 
the SZE-mass scaling relation, and concludes that cooling and 
star formation have a significant effect on the normalization 
but not on the scatter. Moreover, Nelson et al. (2014a} have 
also examined the effect of radiative cooling, star formation 
and energy feedback from SNe on the non-thermal pressure 
fraction and gas velocity anisotropy in group and cluster size 
halos and found little systematic dependence on gas physics. 

To study the evolution of Y, we extract halos from four red- 
shift outputs: z = 0.0,0.5,1.0,1.5. At each redshift we apply 
an additional mass-cut to ensure mass-limited samples at all 
epochs. The mass-cuts and resulting sample sizes are as fol¬ 
lows: 65 clusters with M 200/11 > 6 x 1O 14 /z _i M 0 at z = 0, 48 
clusters with Mioom > 2.5 x 1O 14 /z _1 M 0 at z = 0.5, 42 clusters 
with TLoOm > 1.3 x 1O 14 /z -1 M 0 at z = 1 -0, and 42 clusters with 
Mioom > 7 X 1O 13 /z“>M 0 atz= 1.5. 

3.2. Merger Trees 

To examine the evolution of cluster properties, we track the 
evolution of each cluster’s most massive progenitor. We iden¬ 
tify the main progenitor of each z = 0 cluster by iteratively 
following the dark matter particles at each epoch. Mergers 
are located by identifying objects that cross rsooc of the parent 
cluster. We determine the mass ratio of the merger at an epoch 
when the two systems are sufficiently separated to avoid im¬ 
pacting their respective mass profiles, herein we use the epoch 
at which they initially overlap at a radius of r^ooc- We use the 
ratio of M^ooc at this epoch to identify major mergers, herein 
defined as a mass ratio greater than 1:6. We have also ex¬ 
perimented with other merger ratios, but found that 1:6 gives 
the cleanest trends with the smallest cluster-to-cluster scatter. 
Merging events with mass ratios smaller than this are defined 
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Time since last merger (Gyr) Time since last merger (Gyr) 


Fig. 1.— Mean fractional deviation of thermal F,/, ( red) and F,/, +Y, tl (blue) at (left) and r 2 m c (right) from the best-fit relation for F,/, +F„ f Yf,, , shown as a 
dashed horizontal line in each plot. The shaded regions represent 1 <t deviation from the mean. 


as minor mergers. While this ratio is somewhat lower than 
that used in other studies (e.g. Gottlober et al. |200l) >, careful 
exploration of the merger ratio has shown that this threshold 
best separates mergers that have a strong impact on the mea¬ 
sured ICM properties from those that do not. We examine the 
evolution of ICM properties as a function of time since the 
last major merger, defined as the epoch at which they initially 
overlap at a radius of rsooc- 


4. RESULTS 

Using our simulations, we study the systematic biases in the 
SZE scaling relations due to the changing non-thermal energy 
and cluster dynamical state. We characterize the origin and 
evolution of deviations from the self-similar scaling relation at 
high redshift caused by incomplete thermalization of energy 
from major mergers. 


4.1. Effect of Major Mergers on the Y-M Scaling Relation 

We begin by characterizing the deviation of Y from the best- 
fit relation for K,/, + Y n , at z = 0 (Yf it ) as a function of time since 
the last major merger. In Table [I] we show that Yfi, agrees 
with the theoretically predicted self-similar mass scaling and 
redshift evolution and therefore herein we will use Yf, t as a 
proxy for the self-similar relation. 

Figure IT] shows the fractional deviation of K,/, (red) at rgoo f 
and r" 200 c from Yf,, , plotted against time since the last major 
merger, t mergel , for the 65 z = 0 simulated clusters and their 
most massive progenitors. Each cluster is traced from its lat¬ 
est /merger value at " = 0 back to just before the last major 
merger, and each line in the plot depicts the average Y val¬ 
ues for all clusters as functions of /merger- Note that clusters 
with the same / mer ger in the plot may not be at the same red¬ 
shift. The dashed line depicts Yf,, using the best fit normaliza¬ 
tion for Y t h + Y nt at z = 0 (see Eq. [9), which can be interpreted 
as a line of total energy. During the merger, 0 Gyr < /merger 
< 2 Gyr, we see large deviations from the line of total en¬ 
ergy, a trend also observed by Wik et al. (2008). At this epoch 


the clusters are actively undergoing major mergers, producing 
large merger shocks (see N12 for more details). The shocks, 
which carry a large amount of pressure, result in an elevated 


Y t i,. After 2 Gyr and 3 Gyr at rgoo<? and rjooo respectively, the 
shock has propagated outside of the region and the cluster is 
effectively a closed system. Over the next few Gyr the gas 
motions induced by the merger thermalize and we see a rise 
in the thermal Y by ~ 6% for rgooc and ~ 8% for r 2 ooc- Obser¬ 
vational works on the Y-M relation to date (e.g., Sifon et al. 
2013[ Saliwanchik et al. 2013]) are still too small and noisy to 


reveal this trend, but our results agree with those from simula¬ 
tions by Kay et al. ( |2012| , Battaglia et al. ( 2012[ ) and Krause 
et al. (|2012jl, who observed significant offset between regu¬ 
lar ancfactively merging clusters in their gravitational effects- 
only simulation, where disturbed objects lie slightly below the 
Y — M 500 relation. Specifically, Battaglia et al. ( |2012| ) quan¬ 
tified the dynamical states of their clusters with the ratio of 
kinetic to thermal energy, with a higher value corresponding 
to a more disturbed cluster. These results support our finding 
that there is a significant difference in the fraction of non- 
thermal energy between regular and disturbed clusters, and 
that this difference contributes to the apparent deviation from 
self-similar mass scaling. 

As described in Section |2.3| we calculate the combined 
Y,h + Y nt to estimate the contribution from the missing non- 
thermal energy from the merger-induced gas motions. Y th +Y n , 
is overplotted in blue on Figure [T] For Y measured within 
F 5000 Y,h+Y nt is able to recover the energy conservation in the 
ICM for epochs more than 1.5 Gyr after the merger. As men¬ 
tioned above, the merger shock results in an excess of thermal 
energy that ultimately leaves the enclosed area, resulting in 
an overestimate of the total energy during the merger. For Y 
measured within rgooo there is an excess in Y,/, + Y,„ when the 
cluster is actively merging (/merger < 1 Gyr) and then an under¬ 
estimate in Y,i,+Y nt for the next 2 Gyr. As the shock propagates 
outwards, it carries a significant amount of gas along with it 
to large radius, producing an apparent underestimate before 
the displaced ICM settles back within /inoc This effect can 
also be seen between 1.75 and 2.25 Gyr within rsoor, but to a 
much lesser degree. Beyond 4 Gyr after the beginning of the 
merger, the system returns to a state of energy conservation 
and, by including the non-thermal energy contribution, we are 
able to recover the self-similar mass scaling and evolution to 
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TABLE 1 

Summary of transit parameters 




Y,h 



Y„, + Y n , 


Redshift 

111 A14 

a 

a 

lnA[4 

a 

(7 

0 

-12.01 ±0.07 

1.63 ±0.05 

a 114+0.023 

11 -0.019 

-11 .88 ±0.07 

1.68 ±0.05 

0 115 +0031 

11 -0.023 

0.5 

-12.06 ±0.04 

1.64 ±0.04 

0 103 +0036 
U.iUJ_oo23 

-11.87 ±0.03 

1.69 ±0.03 

0 090 + ° 034 
u.uyu_o oi 8 

1 

-12.12±0.03 

1.71 ±0.05 

0 127+ 0033 
u ’ iZ '- 0.021 

-11.87 ±0.02 

1.76 ±0.04 

0.090t“} 

1.5 

-12.12±0.02 

1.60± 0.05 

0 11 3 +0032 
u - aa:5 -o.022 

-11.84 ±0.02 

1.64 ±0.03 

0.084t°« 

0 (relaxed) 

-12.06 ±0.06 

1.71 ±0.04 

0 070 + ° 007 
u.u /u_o 007 

-11.99 ±0.05 

1.68 ±0.05 

0.058!°;“ 

0 

-12.11 ±0.10 

1.56± 0.05 

0 123+ 0023 

U - Az:> -0.017 

-12.00 ±0.07 

1.62 ±0.04 

0.085!°;“ 

0.5 

-12.28 ±0.06 

1.58 ±0.05 

0 112+ 0027 

11 ^-0.020 

-12.06 ±0.04 

1.63 ±0.03 

0.082!°;°!’ 

1 

-12.43 ±0.04 

1.66± 0.05 

0 1 21 +0 026 
u - 1Zi -0.021 

-12.14 ±0.03 

1.70 ±0.04 

0.085!°;“ 

1.5 

-12.42 ±0.02 

1.59 ±0.05 

0 120 +0 027 

U-A -0.022 

-12.09 ±0.02 

1.61 ±0.04 

o ioo +0 018 

u.iuU-0 016 

0 (relaxed) 

-12.12±0.10 

1.62± 0.06 

n ors+O-Oio 

U.UOJ-Q QQQ 

-12.05 ±0.06 

1.67 ±0.03 

0 Q54+0 006 

u.uj+ o 006 


within 2% for both radii. 


4.2. Effect on the Y—M Scaling Relation 


The above dynamical state dependence and evolution has 
significant implications for the observable-mass, Y — M, rela¬ 
tion. Figure [2] shows the Y-M scaling relations, both with 
and without Y nt , of the simulated clusters at /yoo (left panels) 
and nooc (right panels), with the sample spanning redshifts 0, 
0.5, 1 and 1.5 (shown as circles, triangles, stars and squares, 
respectively). We fit the scaling relation (Eq.[6]» to the samples 
of simulated clusters at each redshift. The dashed line depicts 
Yfj, as fit for the z = 0 sample for each respective Y and radius. 

Table[l]also reports the best-fit slopes and logarithms of the 
best-fit normalizations for scaling relations at each redshift. 
There is no significant redshift evolution in the normaliza¬ 
tion of Y,i, within rjoo, , in agreement with Nagai ( |2006j i, while 
the normalization of Y t h within r 2 ooc shows a deviation from 
the self-similar redshift evolution model towards higher red¬ 
shift. This departure from self-similarity can be clearly seen 
in the green (z = 1.0) and blue (z = 1.5) points, which lie below 
the best-fit line for z = 0. The combination of (a) cluster out¬ 
skirts hav ing larg er fractional non-thermal pressure (e.g. Nel¬ 
son et al. |2014a[ Shi & Komatsu ;2014[ Shi et al. |2014[ ) and 
(b) the expectation that clusters at higher redshift are growing 
more rapidly leads to a large fractional non-thermal pressure 
for a given measured Y at high z, and in turn results in an 
overestimate in mass compared to the self-similar prediction. 
However, when we include the non-thermal energy contribu¬ 
tion from bulk motions, the normalization of K,/, + Y nt is more 
self-similar at all redshifts and radii. 

We also calculate the scatter in Y at a given mass at each 
epoch as 


a = 


N 


S(lny-lnr ,„) 2 

N -1 


( 10 ) 


where N is the size of the sample, Y is either Y,i, or K,/, + Y nl 
as measured directly from the simulation, and Yfi, is the Y 
value predicted by the best-fit relation for each respective Y 
and radius for a given cluster’s mass and redshift, with red¬ 
shift evolution held at the self-similar value. Table [I] reports 
the scatter for both thermal and combined Y. The quoted er¬ 
rors represent the 68% confidence interval, estimated using 
the bootstrap method with 10,000 random samples. The scat¬ 
ter in Y,i, does not evolve with redshift at r 5 oo, , consistent with 
the findings of Sembolini et al. (|2012)> based on their non- 


radiative simulation. We do, however, see a trend towards 
larger scatter in the outer region of clusters, where there is 
a high non-thermal energy fraction due to longer relaxation 
time in the low-density region. At most epochs, the scatter 
in Y r i, + Y„ t , on the other hand, is systematically smaller by 
~ 20-30% than that of Y,i, alone and comparable at r 5f) ( }l . and 

r 200c- 

In Figure [2] we also examine the effects of dynamical state 
on the scaling reaction. Our relaxed subsample (defined using 
clusters at z = 0 that have not experienced any major mergers 
in the previous 4 Gyr) are shown in solid circles. The scat¬ 
ter for this subsample is also given in the last row of Table |T| 
The relaxed subsample always has significantly lower scat¬ 
ter than the entire sample at the same radius and redshift, for 
both Y,i, and Y,i, + Y m . This result agrees with Battaglia et al. 
( |2012| ), who observed that the least disturbed clusters, defined 
as those with the smallest ratios of kinetic to thermal energy, 
show the lowest scatter (11% in their AGN simulation). If we 
compare the scatter of the thermal Y th for the entire sample to 
that of the combined Y,/, + Y„, for the subsample, the redu ction 
in scatter is as large as ~ 50% for both radii. In Section |4.4[ 
we compare the effects of major mergers to more general ac- 
cretion in driving the scatter in the relation. 

4.3. Distribution of Scatter in the Y—M relation 

We examine the relative contributions to the scatter in the 
observable-mass relation from clusters in different stages of 
major mergers. In Figure [3] we show stacked histograms of 
the logarithmic fractional deviations of Y from the self-similar 
predictions for different f mcr „ cr bins, with Y th in the top panels 
and Y t h+Y nt in the bottom panels. The Y values are normalized 
to Yfit calculated as K,/, + Y„, for the z = 0 clusters at each re¬ 
spective radius. The cluster inner regions (rsooc) and outskirts 
(raoorj are shown on the left and right, respectively. The ac¬ 
tively merging systems are shown in red, having experienced 
a major merger in the last 2 Gyr. The most relaxed clusters, 
with ^merger > 4 Gyr, are shown in blue and the intermediate 
stage clusters shown in green. For reference, we have over¬ 
plotted a Gaussian fit, shown as black curves. 

For Y,i, alone, we find that the largest deviations come from 
actively merging clusters (0-2 Gyr), particularly towards 
lower values of Y. Consistent with Figure [T| we find that the 
relaxed subsample tends to lie at the higher end of the scat¬ 
ter distribution at both radii due to the smaller non-thermal 
energy fraction of the ICM and therefore larger Y ,/, value for 
a given mass. Likewise, the intermediate dynamical state is 
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Fig. 2.— The Y — M relations for both Y t f, and Y t h + Y nt at redshifts z = 0, 0.5, 1, 1.5 at rsooc (top) and r 2 ooc (bottom). Top panel of each plot: the relation 
between total mass and the Comptonization parameter, Y. The dashed line shows the best-fit relation at z = 0. Bottom panel of each plot: the fractional deviation 
from the best-fit relation at z = 0. 


between the relaxed and actively merging states. The radial 
dependence of the scatter distribution shown in Tablets also 
clear in top panels, with a wider distribution in the cluster out¬ 
skirts where the fractional missing non-thermal energy contri¬ 
bution is larger. 

When we add Y nt , we see a tightening of the scatter and a 
weakening of the dynamical state dependence at both radii. 
However, the size of the scatter still shows strong dynamical 
state dependence, with the relaxed subsample exhibiting the 
smallest scatter and the actively merging systems the largest. 
There is also a tail in the actively merging distribution towards 
high Y. This is the overestimation shown in Figure |T] due 
to the excess energy from the merger shock during that time 
period. 


4.4. Effect of Other Accretion Processes 

Since clusters do not accrete mass only through major 
mergers, we investigate whether the observed dynamical state 
dependence in the Y — M scaling relation is due primarily to 
major mergers or the general accretion rate, regardless of ac¬ 
cretion mode (including minor mergers and/or smooth accre¬ 
tion). In this section, we generalize our analysis to total mass 
accretion, and eliminate clusters that undergo recent major 


mergers. The remaining clusters allow us to discern the ef¬ 
fect of non-major merger accretion in that period. 

To quantify the total mass accretion of each cluster, we de¬ 
fine the quantity T following Diemer & Kravtsov (2014 1 : 


T = AlogM/Aloga 


( 11 ) 


where a is the expansion factor. Here, T is defined using the 
masses of the main progenitor at z = 0.5 and its descendant at 
Z = 0 within /- 500 c or r 2 oo c - In Figure[4] we show In ( YIYf,,) at 
Z = 0 for each cluster plotted as a function of F for both rsoo,; 
and r 2 ooc- In this figure, Yf„ is calculated separately for F,/, 
and Y t h + Y nt at z = 0 to better highlight the relative distribution 
of scatter. T,;, is shown as red circles, and Y lh + Y nl as blue 
triangles. The clusters at z = 0 without major mergers in the 
last 4 Gyr are highlighted as filled shapes. 

There is a clear trend of Y,u with F, with more slowly accret¬ 
ing systems having higher Y^/Yf,, and more rapidly accreting 
ones having lower Y,h/Yfi,. The lower F clusters are accreting 
mass more slowly than the higher k counterparts and, there¬ 
fore, contain a more thermalized ICM, resulting in a higher 
Y,h for a given mass. This trend is stronger in the outskirts, 
within raoor, where the non-thermal energy fraction in unre¬ 
laxed systems is larger and the clusters at high F have lower 
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Fig. 3.— Histograms of the fractional deviation of Y t h and Y,^ + Y nt from Yf x[ at z = 0 for clusters at r $ooc (left) and r 200 c (right), overplotted with the best-fit 
Gaussian curves. Each histogram is broken down into three ^merger bins: 0-2 Gyr (red), 2-4 Gyr (green) and > 4 Gyr (blue). 
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Fig. 4.— Logarithmic deviation of Y,i, and Y,/, + Y nt at z = 0 from Y fi „ calculated separately for Y t h and Y t h + Y nt , for clusters at rsooc (left) and r 200 c (right), 
plotted as functions of T. The blue triangles and red circles show Y t j x and Y t h + Y nt respectively. Filled symbols are the clusters with ? me rger > 4 Gyr, and open 
symbols represent the remaining clusters. 


values of Y lh . However, when the non-thermal energy contri¬ 
bution is included, Y,i, + Y nt exhibits little to no dependence on 
r at all radii. The vast majority of clusters with large T are 
not highlighted, indicating that clusters that accrete significant 
amounts of mass do so mostly through major mergers. 

There are several clusters with T > 1.5 that are accret¬ 
ing rapidly (representing clusters that have undergone strong 
mass accretion), but no major merger (depicted by the filled 
shapes). However, these clusters lie very close to the self¬ 
similar relation and contribute little to the scatter at large T for 
either Y,/, or Y,h + Y m . The high T clusters with the larger de¬ 
viations from self-similarity are the un-highlighted, recently 
merged systems. Thus, we conclude that the dominant dy¬ 
namical state effects on the scatter in the Y -M relations are 
due to major mergers and other forms of accretion are sub¬ 
dominant. However, it is likely that in the most relaxed sys¬ 
tems, non-thermal pressure injection from ongoing smooth 
accretion and/or minor mergers are responsible for residual 
bias and scatter in the relation. This can be seen in Table |T| 
as the reduction of scatter for the relaxed subsample with the 


addition of the non-thermal energy contribution, especially at 
large radii. 

5. DISCUSSION AND CONCLUSIONS 

In this work, we have studied the relation between clus¬ 
ter dynamical state and scatter and evolution in Y-M scaling 
relation, using a mass-limited sample of massive galaxy clus¬ 
ters from the Omega500 high-resolution hydrodynamic simu¬ 
lation. To characterize the effects of mergers and dynamical 
states of clusters, we analyzed non-radiative simulations in 
which the ICM is heated by gravitational heating processes 
alone and used merger trees to follow the evolution of indi¬ 
vidual clusters during and after mergers and trace their mass 
accretion histories. Our main results are summarized as fol¬ 
lows: 

1. Scatter in the Y-M scaling relation arises largely from 
the evolution of missing non-thermal pressure sourced 
by gas motions in galaxy clusters. The thermal SZE 
signal increases systematically after a major merger, as 
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the energy from random gas motions decays into ther¬ 
mal energy. 

2. We add the non-thermal energy contribution provided 
by random gas motions and find that the scaling rela¬ 
tion for the combined Y exhibits little to no systematic 
evolution and its scatter is reduced by 20-30% at each 
redshift we examined. 


3. The intrinsic log-normal scatter in the Y-M relation 
of a relaxed subsample at z = 0 is about a third smaller 
than that of the entire sample, with a logarithmic l-er 
scatter of only 7.0% and 8.5% compared to 11.4% and 
12.3% in Y,h for the full z = 0 sample at r5ooc and r2ooc 
respectively. 

4. We also examine the effect of other mass accretion pro¬ 
cesses on scatter around the scaling relation, but con¬ 
clude that their contribution is subdominant compared 
to that of major mergers. Thus, we conclude that the 
main effects of dynamical state on scatter are driven by 
major mergers. 


Our results based on cosmological simulations are in gen¬ 
eral agreement with the results of the idealized merger simula¬ 
tions by Wik et al. (j2008| regarding the impact of mergers on 
the SZE signal. However, there are also important quantita¬ 
tive differences that highlight the advantages of cosmological 
simulations. For example, Wik et al. ( |2008| ) reported a steep 
rise (almost by a factor of 2) in Y due to mergers. While we 
also see a similar rise in Figure IT] the increase in Y due to 
merger is more modest with < 80% between the post- and 
pre-merger values. In idealized simulations, a merger exists 
in isolation. In cosmological simulations, a merger occurs 
alongside ongoing smooth accretion and more ubiquitous mi¬ 
nor mergers, which contribute non-negligible non-thermal en¬ 
ergy to the ICM and result in a lower value of Y. 

The post-merger evolution in Y presented herein closely 
echoes the findings in N12 for the hydrostatic mass bias. They 
also find large deviations in the hydrostatic mass estimate due 
to merger shocks for the actively merging systems (/merger < 
2 Gyr) and systematic bias in the mass estimate due to miss¬ 
ing non-thermal pressure for relaxed systems (f mer ger > 4 Gyr). 
Similarly, by accounting for the non-thermal pressure, there 
is a significant improvement in agreement between the hy¬ 
drostatic mass estimate and the true mass. It is interesting, 
although unsurprising, to find such close correlation in the 
timescales related to both the hydrostatic mass bias and the 
scatter in the Y-M scaling relation, highlighting the impact 
and importance of mergers and the unobserved non-thermal 
energy contribution in galaxy clusters. It is worth highlight¬ 
ing that this co-evolution has implications for observational 
studies of the Y-M using the hydrostatic mass estimate. The 
scatter, particularly for unrelaxed clusters, may be affected 
along both axes by the presence of non-thermal pressure con¬ 
tributions. 

There are additional sources of non-thermal pressure sup¬ 
port not included or addressed in our simulation. While the 
inclusion of cooling and star formation does not affect the size 
of the scatter in the Y-M relation (Nagai[2006|, strong energy 
feedback by AGN may enhance the scatter by a few percent, 
especially at high redshift where star formation and AGN are 
active (Battaglia et al. 2012| Fabjan et al. 201 lb . Addition¬ 
ally, magnetic fields and cosmic rays can provide additional 


non-thermal pressure (e.g., Lagana et al.|2010jl, but their con¬ 
tributions are expected to be negligible: the ratio of cosmic 
ray pressure to total pressure is constrained to < 1%, set by 
the 7-ray observations of Fermi -LAT {Fermi -LAT Collabora¬ 
tion |20l4]). Similarly, the typical magnetic field strength of 

< 10/iG in the ICM corresponds to magnetic pressure frac¬ 
tion of < 1%. However, plasma effects can amplify gas tur¬ 
bulence and provide extra non-thermal pressure support (Par¬ 
rish et al. 2012j >, while physical viscosity, on the other hand, 
can lower the non-thermal pressure support by decreasing the 
level of gas turbulence. Moreover, the electron-ion equilibra¬ 
tion timescale is expected to be longer than to the age of the 
universe in the low density galaxy cluster outskirts (r > rzooc), 
which causes the electron temperature and hence the thermal 
SZE signal to be lower than the prediction of hydrodynami- 
cal simulation which assumes the local thermal equilibrium of 
electrons and ions. However, the effect of the non-equilibrium 
electrons on the integrated thermal SZE signal (Y,/,) is of order 

< 3% and < 5% at rsooc and r2ooc> respectively (Rudd & Nagai 
|2009| , Avestruz et al. |2014) >, which is much smaller than the 
effect of merger-induced non-thermal pressure discussed in 
this work. The results presented in this paper based on hydro- 
dynamical simulations serve as a baseline for further studies 
of these effects. 

Recent observational studies of the Y-M scaling relation 
using various observational mass measurements, such as hy¬ 
drostatic mass estimate and gravitational lensing, have found 
an intrinsic log-normal scatter of 20-30% (e.g. Marrone et al. 
|2009| |2012[ Saliwanchik et al. |2013| >. This is twice the size 
of the intrinsic scatter measured from simulations, potentially 
signifying that additional sources of scatter are not accounted 
for in the simulations. While accounting for the non-thermal 
pressure associated with turbulent motions in the ICM signifi¬ 
cantly reduces the scatter, and evolution in the Y-M relation, 
especially in relaxed clusters, gas acceleration can have a non- 
negligible effect, espec ially in unrelaxed systems (Suto et al. 
2013; Lau et al. 120131 Nelson et al. |2014b]>. Focusing on a 
sub-sample of relaxed clusters (which can be defined, for ex¬ 
ample, using morphological measure, such as centroid shift, 
ellipticity, or power ratios, as a proxy for the cluster’s dynam¬ 
ical state) can help minimize the effect of this irreducible bias 
associated with gas acceleration. In practice, there is also ad¬ 
ditional scatter due to projection effects, which could be sig¬ 
nificant especially for low-mass clusters and groups (White 
et al. 2002] Hallman et al. [20071. We will investigate these 
effects in our upcoming work. In observations, one also does 
not know a priori the true values for Ms oo c or nooc as we do in 
simulations. Kravtsov et al. (]2006|) present an iterative algo¬ 
rithm for simultaneously fitting for isnoc and Msooc using scal¬ 
ing relations. Estimating the appropriate radii for the use of 
Y t h and K,/, + Y nt in future observations could also be achieved 
using this technique. Note, however, that any biases in the 
calibration of the Y-M relation can introduce uncertainties in 


Msooc and hence r5oo c . 

Our results have an important implication for cosmological 
constraints derived from SZE cluster surveys. Since Y,i, + Y nl 
is a better proxy for mass than K,/,, our work opens up a possi¬ 
bility to further improve the current robust mass proxy Y,i, by 
measuring and correcting for the missing non-thermal pres¬ 
sure component Ynt. For example, the upcoming ASTRO-H X- 
ray observation satellite will provide the first direct measure¬ 
ment of internal gas motions in clusters, which in turn should 
provide observational constraints on the non-thermal energy 
content in clusters. However, ASTRO-H measurements will 
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be limited to the inner regions (r < = 0.5rsoor) of nearby 

clusters, unless one can dedicate a very long exposure of order 
> 1 Msec (Nagai et al. 2013). Alternatively, the next gener¬ 
ation of high-resolution, multi-wavelength SZE imaging ob¬ 
servations may enable us to measure the non-thermal energy 
component through the kinematic SZ effect signal generated 
by internal gas motions in clusters (Sunyaev et al. |2003| Na¬ 
gai et al. |2003|). Recent theoretical works have found strong 
yet simple relations between the non-thermal pressure and the 
more easily observable thermal pressure profiles, potentially 
supplying an alternate route for estimating the energy contri¬ 
bution from random motions in galaxy clusters (Nelson et al. 
|2014a[ Shi & Komatsu [20l4l Shi et al. |20l4) . 
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